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Peculiar property of electronic order is clarified for the two-channel Kondo lattice. With two conduction electrons per 
site, the order parameter is a composite quantity involving both local and itinerant degrees of freedom. In contrast to the 
ordinary Kondo lattice, a heavy electron band is absent above the transition temperature, but is rapidly formed below 
it. The change of entropy associated with the ordering is found to be close to In \/2 per site. This entropy corresponds 
to the residual entropy in a two-channel Kondo impurity, which has been regarded as due to localized free Majorana 
particles. The present composite order is interpreted as instability of Majorana particles toward non-Kramers conduction 
electrons plus heavy fermions that involve localized electrons. 
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1. Introduction 

Non-Kramers systems with even number of / electrons 
per site as in Pr'^+ and U''+ are expected to show intriguing 
phenomena that is different from ordinary Kramers systems 
with odd number of electrons in Ce- and Yb-based systems. 
Under the cubic crystalline electric field, the localized f^- 
electron state can have a non-Kramers doublet ground state 
where / electrons do not have spin degrees of freedom. This 
doublet interacts through orbital with conduction electrons 
which have spin and orbital degrees of freedom. Thus the two- 
channel Kondo system can be realized.' Experimentally, sev- 
eral candidates such as U2:Thi_a:Ru2Si2 and Prj;Lai_a;Ag2ln 
have been reported.^' ^ Recently the two-channel Kondo sys- 
tems have attracted more attention since the discovery of a se- 
ries of non-Kramers doublet compounds such as PrIr2Zn2o'*'^ 
and PrTi2Al2o.'' 

In the two-channel Kondo impurity, where the localized 
spin is located on a metal with degenerate channel degrees of 
freedom, the scaling theory tells us that the weak and strong 
coupling limits does not become the fixed point. ' The system 
instead chooses the intermediate effective coupling J^s = J* 
as the fixed point, which is in strong contrast with the ordinary 
Kondo model. This scaling flow is illustrated in Fig. 1. The 
ground state has the peculiar residual entropy 5*0 = In >/2 cor- 
responding to a part of localized spin. This remaining degrees 
of freedom is interpreted as the localized Majorana fermions.^ 
In the periodically aligned /-electron systems, on the other 
hand, this partially removed but remaining entropy will lead 
the system to highly non-trivial phase transitions by inter-site 
interactions. 

The simplest description of the periodic non-Kramers dou- 
blets coupled to conduction electrons is given by the two- 
channel Kondo lattice (KL).*^ One of the characteristics in this 
model is the channel symmetry breaking, which is absent in 
the ordinary KL. At one conduction electron per site (quarter 
filling), the possibility of the channel symmetry breaking with 
staggered alignment is pointed out. The stabilization of this 
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phase at finite temperatures is demonstrated in infinite dimen- 
sions. '° At two conduction electrons per site (half filling), on 
the other hand, the dynamical mean-field study has pointed 
out the channel symmetry breaking using exact diagonaliza- 
tion with finite sized bath." 

Recently we have demonstrated that the order parameter 
for the channel-symmetry-broken phase at half filling is the 
composite quantity involving both localized spin and conduc- 
tion electrons.'^ In this paper, we investigate the details of this 
ordered state. By making a comparison with conventional or- 
dered state, we demonstrate that the entropy resolved below 
the transition temperature is related to 5*0 = In \/2 at the non- 
trivial fixed point shown in Fig. 1. We argue that this entropy 
also implies the relevance of the Majorana fermions in the 
two-channel KL. 

The two-channel KL reads 

where Ckaa (cj^acr) annihilation (creation) operator of the 
conduction electron with channel a — 1,2 and pseudo-spin 
(T =t, |. The operators Si and s^a = \ Y.aa' Aaa^<y<y'^ia<y' 
are the localized and conduction pseudo-spins at site i, re- 
spectively. Here cr is the spin- 1/2 Pauli matrix. The Hamilto- 
nian (1) has the S'[/(2) symmetry in both the pseudo-spin and 
channel spaces. The two-channel KL is applied to the non- 
Kramers doublet system; we regard pseudo-spins as orbital 
degrees of freedom, and channels as Kramers doublet origi- 
nating from time-reversal symmetry for conduction electrons. 
Hence we write the symmetry as SU{2)o for pseudo-spin and 
SU{2)t for channels. In the following, we refer to pseudo- 
Jeff = Jeff = J* Jeff = OC 

• > • ^ • 

S'o = ln2 So = lnV2 So = ln2 

Fig. 1. Renormalization flow of the impurity two-channel Kondo model. 
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spin simply as 'spin' except in §6.2. 

For finite-temperature analysis, we have used the dynam- 
ical mean-field theory (DMFT)'^ which becomes the ex- 
act theory in infinite dimensions. The semi-circular density 
of states defined by pa{e) = (2/7r)A/l — (e/Z))^ is taken 
for conduction electrons. We assume the nesting property 
£k+Q — —£k with Q being one of the nesting vectors. 
For numerical method as impurity solver, we employ the 
continuous-time quantum Monte Carlo method (CTQMC).'"^ 
The algorithm we have used is similar to the one explained 
in ref. 15. In the literature, the CTQMC is developed for the 
model where two localized spins interact with two conduction 
bands. When these two localized spins are regarded as identi- 
cal, the model is the same as the two-channel Kondo model. 
Analytic continuation from imaginary Matsubara frequencies 
onto real ones is performed by the Pade approximation. 

The present paper is organized as follows. We investigate in 
§2 the properties of the disordered states. The phase diagram 
of the two-channel KL is shown in §3. Sections §4 and 5 are 
devoted to discussions about the symmetry broken phases. We 
discuss possible applications in §6, and summarize the results 
in §7. 

2. Thermodynamics in Paramagnetic state 

Before we consider the symmetry -broken phases, we study 
basic thermodynamic properties of the paramagnetic state. 
From the transport analysis by the DMFT, Jarrell et al have 
shown that the paramagnetic state is an incoherent metal.** In 
this section, we investigate this phase in terms of thermody- 
namic quantities such as the susceptibility and specific heat. 

2.1 Susceptibilities 

In general, the susceptibility of the operator M. is given by 

X = [{M{t)M^) - {M)^] At, (2) 

JO 

where M.{t) ~ e'^'^Aie^'^^ is the Heisenberg picture with 
imaginary time. As the choice of A^, we consider the charge, 
spin, channel and spin-channel moments of conduction elec- 
trons: 
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where the number of sites is given by N . We have introduced 
the local number operator Uiaa — c^iaa^iaa- We also consider 
the localized-spin moment defined by 

A4pin(g)=iV-'^5fe-'''-«% (7) 

i 

where S'f is the z component of the localized spin. We note 
that / electrons have only the spin susceptibility. 

Figure 2 shows the temperature dependence of the uniform 
susceptibilities x with q = for two cases of the number 
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Fig. 2. (Color online) Uniform spin and channel susceptibilities in the 
paramagnetic phase for (a) nc = 2.0 and (b) nc = 1.5. The insets show 
the local susceptibilities. 
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Fig. 3. Temperature dependence of the specific heat in the paramagnetic 
state at various fillings. The inset is the plot at nc = 2 with linear scale for 
T. More details are discussed in 55.2. 



of conduction electrons = 2 and 1.5 per site. Each x is 
normalized by the respective Curie-like constant C so that 
X — > l/T at high temperatures. At nc = 2 shown in Fig. 2(a), 
the uniform susceptibilities keep constant values even at low 
temperatures. In contrast, local susceptibilities displayed in 
the inset show the logarithmic temperature dependence. On 
the other hand, at nc = 1.5 in Fig. 2(b), the uniform suscepti- 
bilities also show the In T behavior. It seems that the particle- 
hole (p-h) symmetry suppresses the anomalous In T term in 
the uniform susceptibility. The apparent Fermi-liquid like be- 
havior is also observed in the specific heat as discussed next. 
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Fig. 4. (Color online) (a) Phase diagram of the two-channel KL with J = 
0.8 and (b) ordered moments at T = 0.01. In (b), the filled symbols show the 
moment of localized spin, and the blank symbols for conduction electrons, 
(c) Phase diagram near half filling without staggered orders, and (d) channel 
moment at T = 0.01. The full moment in normalized to unity in (b) and (d). 



2.2 Specific heat 

The specific heat C in the paramagnetic state is shown in 
Fig. 3, which is calculated by differentiating the internal en- 
ergy with respect to temperature. The details of numerical 
derivation of the internal energy are given in Appendix. At 
ric = 2, C/T shows a slight peak at T ~ 0.15 and tends 
to a constant for T < 0.06 within the limit of accuracy. Note 
that the fixed-point model (Toulouse limit) of the two-channel 
Kondo impurity shows the same behavior 



C = jT, 



(8) 



where 7 = 7r/(6r) with T corresponding to the resonance 
width of the Majorana fermion.^ The value of 7 is half of an 
ordinary fermion system with the same F. Away from half fill- 
ing, on the other hand, the specific heat shows the logarithmic 
behavior as in the case of uniform susceptibilities shown in 
Fig. 2(b). 

Let us discuss the entropy that remains in the hypotheti- 
cal paramagnetic ground state at half filling. By comparing 
the entropy of an ordered state that has zero entropy at the 
ground state, we can derive the residual entropy. Relegating 
the reasoning to §5.2, we quote the final result for the residual 
entropy 



S{T = 0) ~ 0.55 In 2, 

per site. This value is very close to = 
sponds to the impurity system. 



(9) 

In -\/2 that coiTe- 



3. Phase Diagram and Strong-Coupling Limit 

Now we discuss the ordered states. The phase diagram of 
the two-channel KL for J — 0.8 is shown in Fig. 4(a), which 
is derived from the divergence of the susceptibilities.'^ Here 
we denote the ordering vectors q = and Q by F (ferro) and 
AF (antiferro), respectively. When we take the other values 
of J, there are only quantitative changes. As expected from 
the residual entropy In y/2 in the impurity system, the two- 
channel KL shows ordering in a wide parameter range. Figure 
4(b) shows the order parameters of each phase at T = 0.01, 
which are defined in eqs. (3)-(7). The magnitudes of moments 
behave in a manner similar to the transition temperature. The 




Fig. 5. Schematic illustrations for the ordered states in the two-channel 
KL. Electrons on line show conduction electrons, and the shaded ovals show 
the formation of the Kondo singlet. The filled and open symbols show the 
presence and absence of electrons, respectively. The size of the filled circle 
on the line shows the occupation number of conduction electrons. 



schematic pictures of these ordered states are illustrated in 
Fig. 5. 

The phase around ric = 2 is the antiferro(AF)-spin order 
that has been discussed in ref. 16. We see also the ferro(F)- 
spin ordered phase with q = in the dilute conduction 
electron region. These spin orders are found also in the or- 
dinary KL, and are not specific to the two-channel case. On 
the other hand, we have found the AF-channel order around 
Uc — 1 and F-channel order around ric ^ 1.6. The channel- 
symmetry breaking is a characteristic of the two-channel 
model, and cannot be explained by the RKKY interaction 
since the Hamiltonian (1) includes only the interaction be- 
tween spins. Especially the F-channel order will be discussed 
in greater detail in this paper because of its peculiar and inter- 
esting properties. 

The AF orders in the phase diagram of Fig. 4(a) can be 
qualitatively understood from the atomic limit whose Hamil- 
tonian is given by 



^atomic M ^ ^ ^aa 



JS ■ (Scl + Sc2) 



(10) 



where we omit the site index. The number of conduction 
electrons is adjusted by the chemical potential /i. At ric = 2, 
the ground state is the doubly degenerate states with the total 
spin 1 /2. The degeneracy is lifted by the inter-site interaction 
to cause the AF-spin order On the other hand, the ground state 
at Tic = 1 is the singlet state formed between localized and 
conduction spins. However, conduction spins also have chan- 
nel degrees of freedom, and the ground state is again doubly 
degenerate. The degeneracy is lifted by the inter-site interac- 
tion to cause the AF-channel order. ' At Uc = 0, we have the 
free localized spin at each site. With finite but small concen- 
tration of conduction electrons, they mediates the interaction 
between localized spins to form the F-spin ordering as in the 
ordinary KL. 

We comment on the region around ric ~ 1-5 in Fig. 4(a), 
where the disordered state seems to persist down to T = 0. To 
remove the finite entropy, the system should show some order- 
ing at sufficiently low temperature, but its transition tempera- 
ture can be tiny. Note that the local spin fluctuation remains as 
is clear from the logarithmic temperature dependence of the 
susceptibility in Fig. 2. 
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4. Nature of Composite Order 

4.1 Growth of order parameters 

Interestingly, the picture from the strong coupling limit 
does not work for the F-channel order. In order to investi- 
gate this ordering in detail, we assume that AF orders are 
strongly suppressed. Such a situation is expected in the re- 
alistic system because of e.g. next-nearest-neighbor hopping 
or geometrical frustration. As shown in Fig. 4(c), the transi- 
tion temperature T^y^.^^ of the F-channel phase increases with 
approaching to half filling. The channel moment at T = 0.01 
shown in Fig. 4(d) is finite but tiny. At the maximum transi- 
tion temperature of the F-channel order at ric = 2, the channel 
moment vanishes. Hence, the channel moment is not a proper 
order parameter 

We have identified'^ that a composite quantity defined by 



(11) 



for each site i serves as the proper order parameter In the 
Fourier space, the order parameter is given by ('I'q^o) where 



N' 



(12) 



We note that the composite order is not described by the one- 
body mean fields such as {Si) and {sdi — s^i-i). The expres- 
sion (1 1) is known as the channel-symmetry -breaking pertur- 
bation in the two-channel Kondo impurity.'^ In the present 
case, this term appears by spontaneous symmetry breaking. 

Let us first consider the system at half filling. In addition 
to (5'q=o) the kinetic energy and double occupation become 
dependent on channel below the transition temperature, which 
is first pointed out in ref . 1 1 . In order to characterize these 
quantities, we introduce the normalized order parameter by 



< 1, 



(13) 



where the channel-dependent quantity (a = 1, 2) is cho- 
sen as: 



= < 



(i) 

(ii) 
(iii) 

l(iv) 



(14) 



For the double occupation (iii), one could choose riiafriiai 
instead of eq. (14). However, this expression does not become 
the order parameter under the channel field as will be dis- 
cussed in §4.2. Hence we need to subtract 1/2 from the local 
number operator as in (iii) of eq. (14) for natural definition of 
order parameter 

Figure 6 shows the growth of the order parameter as a 
function of temperature. The composite order parameter (i) is 
the largest among the order parameters, and the kinetic energy 
(ii) and double occupation (iii) have the smaller values. Since 
the coupling between localized and conduction spins is a driv- 
ing force of the order, the primary order parameter should be 
and the resultant imbalance of channels is seen in 
kinetic energies and double occupations as a secondary ef- 
fect. Indeed, as seen in §5.2, the kinetic energy increases by 
the phase transition, which means that the F-channel order is 
disadvantageous for conduction electrons. The interaction en- 
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Fig. 6. (Color online) Normalized order parameter 4> defined in eqs. (13) 
and (14) as a function of temperature. The full moment is unity. 



ergy is lowered by the phase transition. When we focus only 
on the conduction electrons, the quantity (ii) is relevant as the 
order parameter. It is interpreted as the channel ordering in k 
space. 

The order parameter (5'q=o) is a composite quantity, and 
cannot be described by any one-body mean field. Instead the 
effective mean-field Hamiltonian of the F-channel phase is 
given by 



n 



MF 



n 



(15) 



Here A is the magnitude of the two-body mean field. Thus, 
even in the mean-field picture, we still have to solve the many- 
body problem for the F-channel ordering, which is in strong 
contrast to conventional orders. 

4.2 Symmetry analysis 

Respecting the SU (2)t symmetry in the channel space, the 
composite order parameter can also be written in the symmet- 
ric form as'' 
1 



^ ^ ^ ^ Cja^iSli ■ aa')^aa' C-ia' a' — Si ■ ti, (16) 



where we have introduced the spin-channel moment tensor by 



EE"=^"'^'^^o-''^""'C»a'^', (17) 



which can be regarded as a vector either in spin or channel 
space. Especially (t;)^^ has the same form as eq. (6). The z 
component of the vector ^i coincides with the operator de- 
fined in eq. (1 1). It is clear that the order parameter breaks the 
SU{2)t (time-reversal) symmetry. Since SU{2)o symmetry 
is not broken, the orbital Kondo effect can be active in the or- 
dered phase, and can make the pseudo-spin (orbital) singlet. 

In terms of symmetry, the channel moment rrichan, which 
includes mchan as the z-component is also a vector in the 
SU{2)t space. Hence ^ and rrichan are in general mixed. If 
the p-h symmetry is present, however, they are decoupled. Let 
us discuss the situation using the Landau free energy, which 
describes the F-channel order as 



/to, 



(18) 
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Fig. 7. (Color online) Order parameters and channel moment at h = 0.2 
as a function of temperature. The definitions of the symbols are the same as 
Fig. 6. 



where /, B{> 0) and a(> 0) are constants, while A is a linear 
function of temperature. The transition temperature is given 
by the condition A = P/Aa. 

We shall show that / = if the p-h symmetry is present. 
The p-h transformation involving the localized states, is 
defined by 

^Ciaa^ — (^iacn (19) 
^cla^'' ^C^ca, (20) 

^\ia){i<T'\3^^^ = 6aa' - \i(j'){i<j\, (21) 

where \ia) is a localized-spin state at site i. Hence 'J as well 
as the interaction part given by eq. (1) is invariant under J^. 
However we obtain the sign reversal for one-body quantities 
such as the kinetic energy and the channel moment: 

^mchan(g=0)^~' = -mchan(g-O) (22) 

The conduction electrons are p-h symmetric at ric = 2, since 
the chemical potential is at the center of the band. Hence, af- 
ter operating the free energy F remains the same, while 
JTichan changes sign. Then / should be zero in this case. 
Namely, the channel moment may remain zero in the p-h sym- 
metric case. Away from half filling in the F-channel order, the 
tiny channel moment arises as seen in Fig. 4(d). 

The condition I — with the p-h symmetry suggests that 
the spontaneously broken symmetry in the F-channel phase 
may also involve the p-h symmetry. In order to identify the 
relevant symmetry of the composite order, we introduce an 
external field defined by 

Hcxt - -i/i^mcha„(9=0), (23) 

which breaks the channel symmetry, and hence SU{2)t- 
Figure 7 shows the calculated order parameter defined in 
eqs. (13) and (14) for h — 0.2. The channel moment is finite 
even at high temperatures because of the external field. How- 
ever, the composite order parameter becomes finite only 
below Tj,^jan — 0.019. Hence, the symmetry broken below 
the transition temperature is not simply SU{2)t- The rele- 
vant broken symmetry is actually the product where 
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Fig. 8. (Color online) Single-particle spectra A(k,lj) in (a) paramagnetic 
and (b) F-channel ordered states. 



is the time reversal. Note that ^ makes 

^■^g^O-^-' = -*q = 0, (24) 
^mehan (9=0)^-1 = -mchan (<7=0) . (25) 

Since the product leaves mchan invariant, Hcxt in 

eq. (23) does not break the symmetry. This symmetry 
is broken only if 7^ 0, which is odd under 

Away from half filling, generally we have 1^0 since there 
is no p-h symmetry. The symmetry breaking by ^ is simply 
SU {2)t- If we apply Hcxt of eq. (23) away from half filling, 
we have 4* 7^ already at higher temperature because of the 
finite coupling with mchan 7^ 0. Namely the second-order 
transition changes into a crossover. With weak external field, 
however, the channel moment in the ordered state is tiny as 
shown in Fig. 4(d). Hence, the crossover is rapid and should 
look almost like a phase transition. 

4.3 Single-particle spectrum 

Now we consider the single-particle spectrum for the F- 
channel order, which leads us to the effective one-body 
Hamiltonian at low energies. We define the spectrum by 

Ak{uj) ^ --y^lmGkacr{oL! + iri), (26) 

where 77 = +0. The single-particle Green function is defined 
by 

Gkacriz)'^ = Z + ^ - Ek -'>^acriz), (27) 

where z is a complex frequency and Eqct is a self energy. In 
the DMFT, the wave-vector dependence enters to the single- 
particle Green function only through Ek- Hence we introduce 
the parameter k by the relation ~ cos k, and plot the 
spectrum as if the system were in one dimension. 

As shown in Fig. 8(a), the spectrum in the paramagnetic 
state is almost the same as the non-interacting one. The broad 
spectrum is consistent with the incoherent metallic state. ^ 
This non-Fermi liquid behavior originates from interaction 
of localized spins with overscreening degenerate channels. In 
particular the Kondo singlet state cannot be formed. 

On the other hand. Fig. 8(b) shows the result in the 
symmetry-broken phase, where the degenerate condition of 
channels is no longer satisfied. The spectrum near the Fermi 
level is clearly characterized by hybridized behavior of the 
Green function, which has the form 

GkMz) °" (28) 

Z - QaEk - ttaVa /z 
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are also shown. The insets show the magnified picture at low temperatures 
with normal scale. 



Fig. 9. (Color online) Temperature-interaction phase diagram for the half- 
filled two-channel KL. The strong coupling limit J = oo corresponds to 1 in 
the horizontal axis. 



Sa<.(2) = ^+5,z + 0(z2), (29) 

where the renormalization factor is given by = (1 — foa)^^. 
For a = 1, we numerically evaluate ai ~ 0.51 and Vi ~ 0.00 
from analysis of the self energy at T = 0.005. This means 
that the conduction electrons with a — 1 shows the Fermi- 
liquid behavior without hybridization. For a = 2, on the other 
hand, the values are obtained as 02 ~ 0.43, V2 — 0.33. This 
indicates the behavior of the Kondo insulator with effective 
hybridization V2. Thus, the spectrum displays the admixture 
of the Fermi liquid and Kondo insulator. We can understand 
this behavior also from the illustration in Fig. 5(d). In the case 
away from half filling, the Kondo insulator for a = 2 is doped 
into the metallic heavy Fermi liquid at sufficiently low tem- 
peratures. 

The effective model that reproduces the spectrum is de- 
rived from the behaviors of the Green function (28). The low- 
energy behavior of Fig. 8(b) is described by the following 
one-body hybridization model: 

Hcff = ^ ikaci^^Ckaa + V2 ^{ftc^2a + h.c), (30) 
kaa icr 

where ika = o,a£k, and Va = \fa^Va- Here we regard the 
localized spins as "electrons" and introduce the fermion op- 
erator /iCT- From this effective model, it is interpreted that 
the localized spins acquire the itinerancy in the ordered state 
because of the strong renormalization by the Kondo effect. 
Note that this itinerant character is absent in the paramagnetic 
state; we do not see any hybridization of localized electrons 
as shown in Fig. 8(a). This is in contrast to the ordinary KL. 

5. Contrast between F-Channel and AF-Spin Orders 

Peculiar characters of the F-channel order becomes clearer 
by comparing with the AF-spin order at half filling. Roughly 
speaking, the F-channel and AF-spin orders resolve the en- 
tropy In \f2 and In 2, respectively. 



5. 1 Phase diagram 

Let us first discuss the temperature-interaction phase dia- 
gram which includes the information for the origin of order- 
ing. Figure 9(a) shows the T vs J phase diagram for the AF- 
spin order at half filling. In the weak-coupling regime, the 
transition temperature T^^^ scales with as expected from 
the RKKY interaction. For large couplings, on the other hand, 
Tg'^^j is proportional to 1/ J, which means the strong-coupling 
approach account for the order Namely, as discussed in §3, 
the total spin 1/2 states in J —> 00 at each site are split by the 
effective inter-site interaction. Thus the AF-spin order can be 
understood in the above two ways, that is, from either weak 
or strong coupling limit. 

Next we consider the F-channel phase shown in Fig. 9(b). 
In the weak coupling region, the transition temperature scales 
with the Kondo temperature Tk oc exp[— l/po(0) J]. Hence 
the weak-coupling approach fails, since Tk cannot be ex- 
panded by J. For the large coupling region, on the other hand, 
it is likely that the transition temperature rapidly goes to zero 
as J cx), although we cannot go to temperature region 
lower than T < 0.01. This numerical difficulty arises because 
the higher perturbation becomes relevant due to large J. In 
any case, it is clear that the intermediately strong interaction 
is essential for the F-channel ordering. 

The inability to approach the F-channel ordering either 
from weak- or strong-coupling limit is best seen in the single- 
particle spectrum. Namely, in the Fermi liquid channel a = 
1, the effective Kondo coupUng tends to zero, while in the 
Kondo insulator channel a = 2 the coupling tends to infin- 
ity. Thus the F-channel phase is the mixture of weak- and 
strong-coupling Umits depending on channels, and cannot be 
approached from either limit in a perturbation method. This 
feature is a reminiscent of the non-trivial fixed point J* in the 
impurity two-channel Kondo model shown in Fig. 1 . 

5.2 Internal energy, specific heat and entropy 

We compare temperature dependence of the internal energy 
between the two orders at half filling. Figure 10 shows the in- 
ternal energy given by ("H) = (■Hkin) + ('Hint) where 'Hkin and 
Hint are kinetic energy and interaction parts of the Hamilto- 
nian. Note that the internal energy coincides with the expecta- 
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Fig. 11. (Color online) Temperature dependence of specific heat (filled 
symbol) and entropy (blank symbol) for (a) AF-spin and (b) F-channel orders. 
The left and right scales are taken for specific heat and entropy, respectively. 
The specific heat at T = 0.005 in (b) is estimated from the renormalization 
factor. 



tion value of the Hamiltonian because of the condition /i = 
at half fining. 

For the AF-spin order shown in Fig. 10(a), both the kinetic 
and interaction energies are lowered by the phase transition. 
On the other hand, for the F-channel case in Fig. 10(b), ("Hint) 
decreases below the transition temperature, while ('Hkin) in- 
creases. Namely the F-channel symmetry breaking is disad- 
vantageous for conduction electrons. Hence, we confirm that 
the essence of the F-channel ordering is in the interaction 
term. Namely, the quantity (ii) in eq. (14) is not appropriate 
as the primary order parameter 

Next we discuss the specific heat C{T). Figure 1 1(a) shows 
the result for the AF-spin order. The specific heat jumps at the 
second-order transition temperature T^^. Below the transi- 
tion temperature, C (T) goes to zero exponentially, since the 
energy gap is formed at the Fermi level. For F-channel order, 
on the other hand, we have finite values of C{T)/T at low 
temperatures as shown in Fig. 1 1(b). The contribution comes 
from the Fermi-liquid channel a = 1 of conduction electrons. 
The specific heat at low temperatures is estimated from the 
renormahzation factor as 



3ai 



(31) 



Here oi is defined in eqs. (28) and (29). In the present case, we 
obtain 7 ~ 8.2. As indicated in Fig. 11(b), the quasi-particle 
estimate is in excellent agreement with the direct numerical 
evaluation. 

Let us consider the entropy A5 from the ground state up to 
the phase transition. By integrating the specific heat, we de- 
rive the entropy as plotted in Fig. 1 1. In the high-temperature 
hmit, the entropy tends to 5(r 00) = In 32 since the two- 
channel KL per site has 2x4-^ — 32 states with equal weights. 
For the AF-spin order, we estimate as AS" ~ 1.35 In 2 as 
shown in Fig. 1 1(a). This value is interpreted as coming from 
the entropy In 2 of a locahzed spin plus contribution from con- 
duction electrons. 



On the other hand, the F-channel order has AS ~ 0.79 In 2 
as shown in Fig. 1 1(b). The value less than In 2 means that a 
local spin keeps correlation among themselves and with con- 
duction electrons even in the disordered phase. Using this re- 
sult, we now derive the residual entropy in the hypothetical 
ground state. The entropy change in the paramagnetic state 
from r = to T = Tj^^^ is estimated as 0.24 In 2 below 
-^chan' which is illustrated in the inset of Fig. 3(b). Here we 
have assumed that C /T is nearly constant down to the ground 
state. By comparing the entropy just below TJy^^^-^ between the 
hypothetical paramagnetic state and the F-channel state, we 
deduce the residual entropy (0.79 - 0.24) In 2 = 0.55 In 2 as 
given in eq. (9). 

6. Discussion 

6.1 Relation to two-channel Kondo impurity 

As we have seen in §2, the paramagnetic state of the 
two-channel KL with ric — 2 shows remarkable similar- 
ity to the two-channel Kondo impurity system. This can 
partly be understood by the nature of the DMFT. In each 
step of the DMFT iteration, we encounter an effective impu- 
rity model with self-consistently determined local density of 
states. Therefore the convergent paramagnetic solution shares 
the same local property with an impurity system with the cor- 
responding density of states. Note that the specific heat is de- 
termined by the local Green function as seen from eq. (A l). 
Hence it is natural that we obtain the logarithmic temperature 
dependence of the local susceptibility and specific heat for 
ric 7^ 2. What is surprising, however, is that the two-channel 
KL at TT-c = 2 seems to go to the Toulouse limit^ for each 
site that contains free Majorana particles. The linear specific 
heat without logarithmic correction is the clearest evidence, 
and the residual entropy of In \/2 per site strengthens the in- 
terpretation using free Majorana particles. 

We further discuss the relation between the F-channel or- 
der and Majorana particles. According to ref. 7, the relevant 
degrees of freedom for Majorana particles is given by the 
spin-channel moment defined by eq. (17). Thus the compos- 
ite order parameter = Si ■ ti is closely related to Majo- 
rana degrees of freedom. In analogy with the impurity sys- 
tem, the residual entropy in the paramagnetic ground state is 
described by free Majorana fermions centered at each site. 
Then it is tempting to regard the F-channel order as a reso- 
lution of Sq = In \/2 that remains in the two-channel Kondo 
model. Namely, the instability arising from the local Majo- 
rana fermions drives the system to ordering so as to release 
their remaining entropy toward T = 0. 

We have not yet understood why the p-h symmetry brings 
the system toward the Toulouse limit of the two-channel 
Kondo impurity. The origin of nearly constant behavior of the 
homogeneous susceptibility x{T) at low-temperature T may 
also be related to the p-h symmetry. It will be an interesting 
future problem to develop a perturbation theory for the two- 
channel KL from the limit of independent Majorana particles. 

6.2 Channel symmetry breaking described by non-Kramers 
doublet 

We discuss the relevance of our results to the non-Kramers 
doublet systems where the pseudo-spin describes a non- 
Kramers doublet, and the channel index a corresponds to the 
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real spin. Thus the channel- symmetry breaking in the two- 
channel KL implies the breakdown of the time-reversal sym- 
metry. 

Without the p-h symmetry, the F-channel order gives rise to 
a finite magnetic moment corresponding to mchan in Fig- 4(d). 
However, it is too tiny to account for the large thermodynamic 
anomaly in the specific heat. The proper order parameter has 
been identified as the composite quantity which has the 
SU{2)t symmetry but we have considered in this paper the 
^-component. 

The quantity itself does not give rise to a magnetic 
moment, although time -reversal symmetry is broken. Hence, 
is regarded as a magnetic multipole. This multipole mo- 
ment has a radius roughly given by vf/Tk, which is gigantic 
as compared to ordinary multipoles with localized character 
Namely, the order parameter is an itinerant magnetic multi- 
pole centered at each site, and the Kondo effect is responsible 
for the formation. We emphasize that the composite order re- 
quires both itinerant and localized degrees of freedom, and is 
very difficult to treat if one starts from either Umit. 

One may ask whether the composite order is related to the 
hidden order parameter in URu2Si2.^'*'^' The related system 
U2;Thi_j;Ru2Si2 shows the two-channel Kondo behavior,^ 
and URu2Si2 may inherit the property. Then the two-channel 
KL can be a starting point to access URu2Si2. It is likely in 
URu2Si2 that the hidden order is not homogeneous, but has a 
finite ordering vector In order to deal with this situation, we 
need to modify the hopping so as to include the mixing of two 
channels. The higher-rank magnetic multipole in the present 
paper should be hard to be detected by any diffraction exper- 
iment although the system shows the large thermodynamic 
and transport anomalies. Hence the composite order param- 
eter is a candidate of the hidden order parameter Under 
the fourfold crystal symmetry, the transverse components 
or form a two-dimensional Eu representation. The dou- 
bly degenerate representation seems consistent with recently 
performed experiments.^^' 

6.3 Formation of narrow bands by ordering 

As seen in Fig. 8, the narrow bands are newly formed by 
the phase transition into the F-channel ordered state, which 
is not seen in the paramagnetic phase. Such behavior is ob- 
served also in the other ordered states such as AF-spin and 
AF-channel phases in the two-channel KL. 

Let us compare this characteristic temperature dependence 
of the spectrum with that in the ordinary KL.^"* In the ferro- 
magnetic KL, where the Kondo effect does not operate, the 
clear narrow band is never observed. The antiferromagnetic 
KL, in which the Kondo effect is working, heavy bands are 
formed in the ordered state. However, the clear hybridized 
feature of electronic spectra has already been observed in 
the disordered state. This is the principal difference between 
single- and two-channel KLs. This difference originates from 
the difficulty of forming a Kondo singlet in the two channel 
model. Only by ordering, the localized spin can find a proper 
partner among the two. 

The formation of narrow bands by ordering is recently ob- 
served in URu2Si2 by the single-particle spectroscopy. In the 
angle-resolved photoemission spectroscopy (ARPES),^^'^* 
the flat band appears below the transition temperature. On the 



contrary, it is not the case for the antiferromagnetic transition 
in a Rh-substituted sample of URu2Si2. Hence the appearance 
of the electronic band is a characteristic of the hidden order. 
A similar behavior is also observed in the scanning tunnehng 
spectroscopy (STS).^^"^^ We postulate that this formation of 
narrow bands by ordering can be explained by the mechanism 
presented in this paper. 

7. Summary 

With use of the DMFTh-CTQMC approach, we have in- 
vestigated the two-channel KL as the simplest model that de- 
scribes periodic non-Kramers doublets coupled with conduc- 
tion electrons. In the paramagnetic state, the homogeneous 
susceptibiUty x and specific heat C at half filUng show the ap- 
parent Fermi liquid behavior Away from half filling, without 
the p-h symmetry, both x and C show logarithmic tempera- 
ture dependence. We have derived the entropy at half fiUing, 
and found that nearly In \/2 per site remains in the paramag- 
netic ground state. 

The channel symmetry breaking at half filUng is pecu- 
liar in that the ordered state is characterized by the prod- 
uct of p-h symmetry !^ and the time-reversal symmetry ^. 
The channel-asymmetric composite quantity between local- 
ized and conduction spins changes sign under this product, 
while the channel moment remains intact. Thus the broken 
symmetry is identified as the product The electronic 

spectrum shows the drastic change at the transition temper- 
ature, where the narrow bands are formed only in the or- 
dered phase. By comparing thermodynamic quantities in the 
F-channel and the AF-spin phases, we have proposed that the 
F-channel and AF-spin orders are characterized as resolution 
of entropy In \/2 and In 2, respectively. 

In analogy with the two-channel Kondo impurity, the re- 
maining entropy = In in the paramagnetic state is de- 
scribed by the localized Majorana fermions. The composite 
order is regarded as the change of these Majorana fermions 
into ordinary fermions by phase transition, since the ordered 
state is naturally described by the Fermi liquid theory. In this 
paper, we have considered the Hamiltonian (1) with the sim- 
ple isotropic Kondo exchange coupling. However, the form 
of interaction becomes more complex in general cases. For 
example, the Fs doublet in the tetragonal symmetry as in 
URu2Si2 has a dipole moment in addition to quadrupole mo- 
ment. Another probable situation URu2Si2 is that the two sin- 
glet levels are almost degenerate and form a pseudo doublet. 
It is an interesting future work to study how the behaviors 
obtained in the two-channel KL are modified by these com- 
plications. 
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Appendix: Numerical Evaluation of Internal Energy 

We calculate the expectation value of the Hamiltonian from 
the local Green function as 



(:W)=TVi£„Gehg(ien)e'^''^ 



(A-1) 



where e„ = (2n + l)7rT is the fermionic Matsubara fre- 
quency. We have defined the "charge" Green function by 
Gchg = J2aa ^oia where 



is the local Green function. The internal energy E is given by 
the relation {V) = E - iiN. 

Owing to the presence of the form factor ie„ in eq. (A l), 
we have to numerically sum up a large number of Matsubara 
frequencies. In order to reduce this difficulty, let us consider 
the expansion from high-frequency Umit 

A B 



Gchg(i£n) — 



AGchg(i£n), (A-3) 



(ie„)2 (ie„)3 

where c — Y^^^ 1 = 4. The coefficients A and B are real. 
The function AG is introduced so that it behaves as (i£„)~^ 
when £„ — >^ oo. Then we obtain the relation 



(^) = ^ - ^ +T5^i£„AGehg(i£„). 



(A-4) 



Note that we do not need the convergence factor in this 
expression. In the numerical calculation, the summation in 
eq. (A-4) is performed in the range \n\ < L. We typically 
take L ^ 10^ in the simulation. The real coefficients A and B 
are numerically evaluated from the high-frequency behavior 
of the Green function as 



A= lim Re(ie„)^Gchg(ien), 



B = lim Re (ie„)^ 



Gchg(ien) 



(A-5) 



(A-6) 



We have to carefully determine A so that AGchg does not 
include the (ie„)~^ term. Otherwise, numerical errors can 
arise, since we need a convergence factor in the Matsubara- 
frequency summation in eq. (A-4). Especially in the p-h sym- 
metric case, Gchg(i£n) is purely imaginary, and hence the co- 
efficient A is exactly zero. 
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